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Abstract 

We discuss the Hartree equation arising in the mean-field limit of large systems of bosons 
and explain its importance within the class of nonlinear Schrodinger equations. Of special 
interest to us is the Hartree equation with focusing nonlinearity (attractive two-body interac- 
tions). Rigorous results for the Hartree equation are presented concerning: 1) its derivation 
from the quantum theory of large systems of bosons, 2) existence and stability of Hartree 
solitons, and 3) its point-particle (Newtonian) limit. Some open problems are described. 



*Notes to a lecture presented by J. F. in the PDE seminar at X-Palaiseau, Paris, April 2004. 
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1 Introduction 



We present an overview of rigorous results and open problems concerning the Hartree equation 



where — ip(t, x) is complex- valued, and V and are real- valued functions on R"^ . Furthermore, 
A denotes the Laplacian and * convolution of functions on M.^. 

Equation (|HEp has many interesting applications in the quantum theory of large systems of 
non-relativistic bosonic atoms and molecules. In fact, this equation arises in the study of the mean- 
field limit of such systems, i.e., of a regime where the number of bosons is very large, but the 
interactions between them are weak. It is therefore of considerable interest to extend mathematical 
methods originally developed for nonlinear Schrodinger equations with local nonlinearities to the 
study of equation (|HE|I and to develop new methods. 

More specifically, the work reviewed in these notes has been motivated by the following physical 
questions. 

(A) What is the behavior of many bosons with attractive two-body interactions in the mean-field 



(B) How does quantum mechanics reproduce Newtonian mechanics of point-particles in limiting 



In these notes we review results providing partial answers to these questions. Our discussion 
is arranged as follows. 

• In Section |21 we discuss the mean-field limit of a bosonic system and the emergence of 
the Hartree equation (|HE|I in this limiting regime. We present rigorous results, physical 
applications, and two open problems related to this limit. 

• Section 13 contains a short review of the Hamiltonian structure connected with (|HE|I . We 
also review some conservation laws and various identities. 

• In Sectional we study solitary wave solutions of equation (|HE|I . which we call "Hartree 
solitons". The main mathematical results are summarized in Theorems |21 and |3| asserting 
existence and orbital stability of Hartree ground state solitons for attractive, "subcritical" 
potentials $. These results are proven by using variational methods. Uniqueness of Hartree 
ground states in general situations is still open and posed as a problem. At the end of this 
section, we address the delicate issue of asymptotic stability and suggest an open problem. 

• In Section[3 we discuss the point-particle (Newtonian) limit of l|HE() . That is, by considering 
a suitable limiting process, we derive Newton's equations of motion for a system of point-like 
particles corrected by a term representing "radiation damping" . 

• In Sectional we draw attention to some further open problems, and offer some conclusions 
and an outlook. 

Notation and Units. In what follows, physical space is chosen to be R^, but the discussion 
can be extended to M'', for d > 3. Throughout this text, Lebesgue spaces of complex-valued 
functions on M.^ are denoted by iP(R^), with norm ||-|lp, 1 < p < oo. The Sobolev space 





(HE) 



limit? 



regimes? 
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H^{M.^) := {ip e L^{R^) : ||VV-'||2 < 00} plays an important role. Physical units are chosen such 
that Planck's constant h = 1. 
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2 Mean-Field Limit 

We consider a system of N spinless, identical, non-relativistic bosons with two-body interactions 
given by a potential, k^{x), where k > is a coupling constant. These bosons are also subject 
to a time-independent, external potential, XV (x), with coupling constant A > 0. If the mutual 
interactions among these bosons are weak, in a sense to be made precise below, the system can 
be interpreted as a quantum Base gas under the influence of an external potential called "trap" 
in condensation experiments. 

For simplicity, we only consider pure states of this A^-particle system. The dynamical evolution 
of a pure state, n, is governed by the linear Schrodinger equation 

idt^^S^N = -ffjv^-jv , (2.1) 
where — ^7v(^, xi, . . . , with Xj £ M.^, ^ < j < N. The A^-particle Hamiltonian is given by 

N N 
j = l i<j 

where (g)s denotes the symmetric tensor product. In equation H2.2|) . denotes the Laplacian 
associated with the jth copy of R^; the boson mass is m. 

Since the total energy of the system scales like 0{N) + kO{N'^), the energy per particle is 0(1) 
if the coupling obeys k — 0{N^^). Thus, one defines the mean-field limit to correspond to 



N ^ 00 and k ^ , such that v :— k. ■ N — const. 



(2.3) 



As shown below, this limit leads to the Hartree equation as an effective description of the dynamics 
of so-called coherent states, or condensates, of the system. 

Remarks. 1) Note that the statistics of bosons dictates this scaling, whereas a fermionic system 
would require a different scaling k = 0{N~^^^). 

2) There is a second limiting regime referred to as "classical limit". It is constructed as follows. 
Approach the mean-field limit, rescale time according to t — y^r, and take m — > 00; (this is 
equivalent to ft ^ 0). In this limit, the Vlasov equation emerges from (jniES; see |(;iVIPn3| and 
references given there. 



2.1 Heuristic Discussion 

We now demonstrate heuristically the emergence of the Hartree equation in the description of 
a bosonic system in its mean-field limit. For this purpose, we investigate the limiting regime of 
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a condensate where all - except for o(N) - bosons arc in the same one-particle state, which is 
described by a wave function tjj E L^(R'^). We choose a coherent state as an initial datum for H2.1|l . 
That is, we pick -00 € L-^(R'^), with HV'olli — li introduce the iV-particle state 

N 

*jv(i = 0,xi,...,xn) = YIMx]). (2.4) 

i=i 

which apparently belongs to i^(R'^)®=^. Assuming that the linear Schrodinger equation (|2.1(l 
approximately preserves the coherent nature of ^'Ar(i) when N becomes large, we can write 

N 

)^l[ij;{t,x,). (2.5) 

Physically speaking, this means that correlation effects remain small. 

Approaching the mean- field limit, we expect that the potential per particle is given in terms 
of an effective mean- field potential, Ueff , as follows 

v^e{t, x) := XV{x) + ^ V / ^{x - Xj)\ip{t, Xj)\^ dxj 

^ ^iJv' (2.6) 

= XV{x)^v{<^*\ilj{t)\'^){x), 

where 

{f*9){x) / f{x-y)g{y)Ay. (2.7) 



From this heuristic discussion we conclude that the dynamical evolution of the bosonic system in 
its mean-field regime is described by the Schrodinger equation for the one-particle wave function, 
ip{t,x), with a potential term given by Vcff(t,x). Thus, we are led to the (nonlinear) Hartree 
equation as an effective description of the limiting dynamics: 



idt%l) = -Alp + A^V + !^($ * IV'HV' , 
ipiO,x) = il^o{x) , 



(HE) 



and we will choose units such that 

TO ==1/2, (2.8) 

from now on. In l|HE|) . the complex-valued function, ip — 'ilj{t,x), with x E M^, describes the 
one-particle state of the condensate with initial state "0(0,2:) = "00 (2;)- For convenience, we keep 
the coupling constants, A > and v > 0, explicit throughout the following. 



2.2 Rigorous Results 

A derivation of the Hartree equation in the mean-field limit of many bosons has been achieved 
with mathematical rigor. The most delicate step to be taken in this analysis is to show that 
correlation effects actually do remain small for all times, i. e., approximation (|2.5|l is justified. 

The first rigorous derivation of l)HE(l dates back to K. Hepp's work |Hep74| , where $ is assumed 
to be smooth and bounded. A generalization valid for all bounded potentials and one-particle 
density matrices (thus admitting mixed states) is due to H. Spohn |Spo80| . Recent work by 
L. Erdos and H.-T. Yau |EY01| covers a rigorous derivation for Coulomb potentials <i>(x) = ±|a;|~^. 
New simple proofs will appear in |Sch04| . 
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We briefly sketch some of the main ideas in the rigorous works cited above. For this purpose, we 
introduce m-particle operators, Am, defined in terms of integral kernels, a(xi, . . . , Xm] yi, ■ ■ ■ , Vm), 
which are assumed to be smooth, bounded functions on M^'" x R'^™, 1 < to < iV. The expectation 
value of Am in an A^-particle state, ^'tv, is given by 



(*Ar,^„*Ar) = / M {xi, . . . , X n) a{xi, . . . , Xm] Vl, ■ ■ ■ ^Vm) 

JV (2.9) 
^'jv(yi, . . . ,yAr) Jl 5{xj ~yj)dxdy . 

j—m+l 

If the A'^-particle initial state is coherent, i.e., 

N 

^N{t = 0, xi, . . . ,xjv) = n ^o(^j) ' (2.10) 
for some ■(/'o G L^(K'^), then, under certain assumptions on V and V'O) 

lim {^N{t).,Am'^N{t)) = I a{xi, . . . ,Xm;yi, ■ ■ ■ ,ym) 



Y[tk{t,Xj)tlj{t,yj)dx dy, 



(2.11) 



as proven in the papers quoted above. Here ^'Ar(t) denotes the solution of the linear Schrodinger 
equation l|2.1(l with initial condition H2.10|l . whereas tp = 'ip{t, x) solves the Hartree equation (|HE|I 
with initial datum tpQ. Equation H2.11|l justifies our heuristic discussion of Section lTTl All relevant 
information on the dynamics of coherent states of the bosonic system in its mean-field limit can 
be obtained from the Hartree equation. 

For our mathematical analysis of the Hartree equation it is convenient to introduce a class of 
potentials defined as follows. A potential $ is called subcritical if 

$ e L'P{R^) + L°°{R^) , for some p > 3/2 , (2.12) 

and likewise for V. Furthermore, we define the energy space for IjHEjl to be 

X {^j e H\R^) : V\ij\^ e L\R'^)} , (2.13) 

equipped with the norm \\4>\\x := HV'b + llV^Alb + || iV^I^/^V^Ha. 

Referring to |Caz03j and references therein, we can state the following result concerning the 
Cauchy problem for IjHEfl in energy space. 



Theorem 1. [Global Well-Posedness] Let $ be real-valued, radial, and subcritical. Assume 
that V is real-valued and satisfies one of the following assumptions. 

a) V is subcritical (i. e., V is a "localized trap", orV = 0). 

b) V € C°°{R^), V >0, and \D°'V\ G L°°{R^) for all \a\ > 2 (i. e., V is a smooth "confining 
trap" with at most quadratic growth). 



Then \HE^ is globally well-posed in X. That is, for every initial datum tpo G X , there exists a 
unique solution tp € C{R,X) depending continuously on ipQ. 
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Remarks. 1) Note that if condition a) holds, then we have that X = H^{M.'^) by virtue of Sobolev's 
inequahties. 

2) The study of the Cauchy problem for (|HE|I dates back to the work of J. Ginibre and G. Velo 

fnvsoi . 

3) We expect that global well-posedness for appropriately chosen initial conditions can also be 
proven for external potentials growing more rapidly than as oo; cf. Theorem El below. 

4) The conditions of Theorem ^ also imply that the A^-particle Hamiltonian, Hn, is bounded 
from below due to the Kato-Sobolev inequalities. In addition, Hj^ then defines a unique selfadjoint 
operator via its form sum; see, e. g., 

In the following discussion, we mainly consider attractive, subcritical ^'s. We derive existence 
of Hartree ground states and establish their orbital stability for such two-body potentials; see 
Theorems 121 and 13 below. We will, however, also comment on the critical case, ^(x) = — 
which misses to be subcritical and leads to blow-up phenomena in IjHEj) of interest in physics. 

2.3 Physical Applications 

Here we describe three physical examples where the mean-field limit yields an appropriate descrip- 
tion. 

Repulsive Interactions 

Bose-Einstein condensation (BEG) in gases with very weak repulsive two-body interactions can be 
found in systems of ^^Rb or ^^Na atoms. In this repulsive case, the approximation of $(a;) by the 
delta function, +5{x), is adequate, and (|HE|I becomes a cubic nonlinear Schrodinger equation of 
local type, i. e., 

idt^p = -Aij + V^p+\ip\^ij, (2.14) 

which is often referred to as the Gross-Pitaevskii equation; see |LSSY03l for a comprehensive re- 
view of rigorous results on BEG with repulsive interactions. 

Attractive Interactions 

BEG in gases with very weak attractive two-body interactions is observed, for example, in systems 
of ^Li atoms, as long as the gas in the trap has a sufficiently low density. But even for A = there 
is a critical number, > 0, such that if 

kN > (2.15) 

bound clusters of atoms exist. This feature is reflected in the existence of Hartree solitons when 
no external potential is present; see Theorem|2 below, for a precise statement. 

In experiments, a collapse process is observed if the number of ^Li atoms in the trap exceeds 
some threshold value. But we remark that, in a theoretical description of such collapse processes, it 
is no longer justified to treat bosonic atoms as point-particles. Their internal structure, ionization 
processes and recombination, and interactions with the electromagnetic field are important when 
the density of the system becomes large. Neither the Gross-Pitaevskii equation (with focusing 
nonlincarity) nor the Hartree equation take such processes into account. One may expect that, at 
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best, the Gross-Pitaveskii equation (with focusing nonhnearity) and the Hartree equation (with 
$ = an attractive van der Waals potential) give a quahtatively satisfactory account of the onset 
of the coUapse process for kN > v^,. 

Boson Stars 

Boson stars are objects, at least theoretical ones, made of bosons subject to their own gravity. It 
is reasonable to use mean-field theory to describe such objects. A model of a boson star in its 
non-rclativistic regime is given by 

i9tV = -AV-- i^da^r^ (2.16) 

which is also referred to as the Schrddinger-Newton equation (or the Choquard equation in another 
physical context; see |Lie77| '). Existence and stability of Hartree solitons arising from ground 
states can be shown for H2.16|l : see Theorems |21 and |31 below. We refer to |HMT03| for a numerical 
study of ifTT?^ . 

An improved model of a boson star arises when high velocities of the bosons are taken into 
account by incorporating special relativistic effects, but still assuming Newtonian gravity. That 
is, we consider the semi-relativistic Hartree equation 

idti) = (V-A + m2 - m)i} - v{\x\-^ * IV-HV' , (2.17) 

with H'^/'^i^) as energy space. Equation H2.17|l leads to a Chandrasekhar type theory of boson 
stars. See |LY87| for a study of the time-independent equation; an extended discussion of H2.17|l will 
appear in |Lenn4| . Note that, since the relativistic kinetic energy has a dimension of length^ ^, the 
two-body potential, <i>(x) = — is critical, large initial data can lead to blow-up phenomena, 

which can be thought of as an indication of gravitational collapse. Of course, this model has its 
limitations. Nevertheless, it is useful to explore the onset of collapse and structure formation of 
bosonic matter by studying equations such as (|2.17|l . 

2.4 Two Open Problems about the Mean-Field Limit 
Mean-Field Spectral Analysis 

As mentioned above, the Hartree equation captures the dynamics of a bosonic system in the mean- 
field limit. We now ask how the bound state energies of the A^-particle Hamiltonian, Hm, can be 
derived from an effective one-particle Schrodinger operator, as the mean-field limit is approached. 

To gain some insight into this problem, we propose the following simplified "effective" descrip- 
tion. We define the ground state energy of Hjs; by 

E% ■.= inia{HN). (2.18) 

This quantity can be related to the Hartree equation ljHE|) as follows. On the basis of results in 
|LY87| . we expect that 

E% = Neo + o{N), (2.19) 
as — > oo, with v ^ n ■ N fixed, where 

eo := mi{2H[ij] : V e A, ||^||2 = l} , (2.20) 
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and 7i [■(/'] denotes the Hartree Hamilton functional; see H3.2|l below. 

As shown in Theorem |21 below, there exists a minimizer, ^ — Q, for problem (|2.2Q(I . which is 
called a Hartree ground state and satisfies the Euler-Lagrange equation 

-l^Q + \VQ + v{^*\Q\^)Q^eQQ. (2.21) 

Next, we propose to generalize equations H2.19I I2.2Q|I to a relation between the bound state 
energies of Hm and the spectrum of an effective one-particle Hamiltonian given by 

iJ^ff := -A + AT/ + z/($ * IQP) • (2.22) 

For simplicity, we assume that Hcs and Hjsi have finitely many bound states (after splitting off 
the center-of-mass motion of Hm if A = 0). We consider the eigenvalue problem 

-ffcff'/'i = Ei^i , (2.23) 

with {eo, ei, £2, • . .} the eigenvalues of iJeff and {4>q = Q,4>i^4'2, ■ ■ ■} the corresponding eigenvec- 
tors. The mean-field picture leads to the following conjecture: Every bound state energy of Hn 
below its contimmm threshold, Sjv, is given by the ground state energy of H^-j plus a suitable 
combination of j energies from the set of e^'s. More precisely, we expect that the discrete spectrum 
of Hm is given by 

ad{HN) ^i^E -.E^ E%_^ + ^(e,, - £0 + eo) + o(l) < Sat , < j < ivj , (2.24) 

with Tin = -E'at-I' as — > 00, with v = k ■ N fixed. Thus, we propose the following problem. 

Problem 1. Show that, for a class of two-body potentials \2.24\! holds as the mean-field limit 
is approached. 

RelELxation into Coherent States 

Next, we describe the problem of "relaxation into coherent states" motivated by our expectation 
that coherent bosonic condensates emerge from non-coherent initial states. 

Recall from equation H2.10|l that derivations of (jHE|l require that the initial A^-particle wave 
function, ^'7v(0), is coherent. Obviously, such initial conditions are exceptional. Therefore it is of 
fundamental interest to understand how a state, ^'Ar(t), arising from an arbitrary initial condition, 
^Ar(O), approaches (in a sense to be made precise) a coherent state in the limit 

t ^ 00 and N ^ 00 , v = k ■ N — const. , (2.25) 

where t will depend on N . Thus, we are led to the following problem. 

Problem 2. Show that, for a class of two-body potentials and a suitable limiting process l^2.25\l . 
non-coherent states relax into coherent states. 

Remark. The study of Problem|21necessitates a mean-field spectral analysis of Hn that also takes 
into account resonances and their life-times. 
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We sketch two different situations, depending on the external potential, and their conjectural 
qualitative behavior in the limit (|2.25() . 

Trapping. For localized traps, XV, physical intuition suggests a relaxation process correspond- 
ing to the "trapping of a Hartree soliton". The TV-particle state, ^'Ar(f), is expected to approach 
(with positive probability g > 0) a, coherent bound cluster of bosons described by a Hartree ground 
state soliton, plus some dispersive "radiation", in the limit l|2.25|l . See Figure Q for a sketch of 
this situation. To arrive at a formal statement for this relaxation process, we use the notation of 
Section We expect that 

f " 
(*jv(0:-4™*Ar(i)) ~ g a{xi,. . . ,Xjn;yi, ■ ■ ■ ,y7n) \\Q{xj)Q{yj)dxdy , (2.26) 

where < g < 1, and Q is a Hartree ground state with \\Q\\2 ~ 1, and where is replaced by 
> in (|HE|I . with iVi < N. See Section^ for details on Hartree ground states. 




Figure 1: "Trapping of a Hartree Soliton". The A^-particle state, '^^{t), approaches a coherent 
state corresponding to a Hartree ground state, Q, which is located near the minimum of the trap. 
Dispersive "radiation" is emitted during this relaxation process. 

In addition to such trapping processes, another situation can occur: If ^'Ar(t) "scatters" off the 
trap, then the A^-particle state can exhibit the same behavior as described next. 

Vanishing external potential. Assume that A = 0, i.e., no external potential is present. If the 
two-body potential is sufficiently attractive (i.e., > i'*), then we expect 4'Ar(i) to evolve into 
bound coherent clusters of bosons receding from each other, plus some small dispersive "radiation" , 
as t ^ cx); see Figure El Using the notation of Section the situation envisaged here can be 
described as follows. Given an m-particle operator. Ami we define Am{Y) via its integral kernel 

a{xi^Y,...,Xrn-Y;yi-Y,...,ym-Y), for F G _ (2.27) 

In view of ()2.11|l . we expect that 

/ (*Ar(0,A„(n*ivW>dr Va- V [ dY [ a{xi-Y,...,x,n-Y; 

JR3 . - , JR3 Jr3™xR3" 

(2.28) 

2/1 - y, ... , vrn -Y)fl {t, x.Ovf {t, yj) dx dy . 

Here il)^^\t) solves HHE|) . with ||V'fc^(0ll2 = 1 and Vi — j^v > v^, and gk are non- negative weights 
with "Y^f, Qk ^ ^- By physical reasoning, we expect "ground state selection", i.e., the one-particle 
wave functions, ip^^^t), asymptotically approach Hartree ground states. 
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t — > oo 



*jv(i = 0) 




Figure 2: An TV-particlc state, '^^{t), evolving, as i — *■ oo, into fc > 1 bound coherent clusters 
of bosons with > 
dispersive "radiation" . 



of bosons with > v^,, i — 1, . . . , fc, receding from each other with velocities, Vi, plus some 



3 Hamiltonian Structure 

Rephrasing partial differential equations as Hamiltonian equations of motion on some infinite- 
dimensional phase space is a well-known procedure. For the Hartree equation, this also clarifies 
its point-particle (Newtonian) limit, which is presented in Section |31 below. 



3.1 Basic Setup 

We identify the phase space for 1HE|) with its energy space, i. e., we put 

r := X = {tA e -ff^(M3) : e L^{R^)} , 



(3.1) 



where it is convenient to treat and ^l> as canonically conjugate, independent variables. Defining 
the Hartree Hamilton functional by 



^ ./B3 L Z 



dx . 



we are led to the Hartree equation as a Hamiltonian equation of motion: 



(3.2) 



(3.3) 



where 6 /Sip denotes the Frechet derivative with respect to ip. Alternatively, (|3.3fl can also be 
expressed in Liouville form 

m = [nim.mw)]^^, (3.4) 

with Poisson brackets defined by 



{V'(x), i/'(y)}pB = {?/'(a;), i/'(y)}pB = , {ipix), V'(2/)}pB = 2i(5(a; - y) . 



(3.5) 
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We remark that if V and $ meet the conditions of Theorem ^ then it is straightforward to show 
(using Sobolev's inequahties) that 7i is a C^-functional on T bounded from below on the spheres 

Sat := {Ver:AA[V;,V] =iV} , (3.6) 

for aU iV < oo, where we define the charge (or mass) functional by 



(3.7) 



Hamilton's equations of motion (|3.3|1 can also be derived from an action principle. The corre- 
sponding action functional is given by 



dt, (3.8) 



and is defined for all C^-curves [ti,<2] ^{t) <= T. The Hartree equation is equivalent to the 
Euler-Lagrange equation for 5[V','0] varied with respect to tp, i.e., 

^ Sm-),i;{-)]=0, (3.9) 



subject to the boundary conditions 

S^Pit)\t=t, = , with J = 1,2. (3.10) 

Note that the canonically conjugate variables, V' and V'l appear in the action principle, i.e., we 
consider a Hamiltonian action principle on phase space. 

3.2 Conservation Laws and Identities 

From now on, we abbreviate our notation and write 

n[ij]=H[^P,^], AA[^]=AA[^,V], S[^P] = S[^i-),^P{-)]. (3.11) 
Charge and Energy 

The invariance of the action. Slip], with respect to U{1) gauge-transformations and time transla- 
tions leads to the conservation of charge and energy, respectively. We have that 

u[m] = , mm] - n^Po] , (3.12) 

where ip{t) solves (|HE|I . In the setting of Theorem these conservation laws hold rigorously and 
play a central role in its proof. 

Linear and Angular Momentum 

If XV = then further conservation laws derived from the Galilean invariance of the action S[ip] 
arise. Because of translation invariance we have conservation of linear momentum: 

V[tp{t)] ^ V[2po] , where V[tp] := -i [ ^\/ipdx. (3.13) 
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Moreover, the center- of-mass motion is inertial, which is expressed by the fact 

■4}{t){x + 2itV)->p{t)dx = const. (3.14) 



Note that, for H3.14|l to hold rigorously, one has to require suitable decay of the initial datum, as 
|a:| oo. 

Finally, from the invariance of with respect to spatial rotations we can infer conservation 
of angular momentum: 

C[ip{t)] = C[i;o] , where C[ip] := -i / ipix x \')^dx. (3.15) 



Various Identities 

We note that there are useful "identities" for nonlinear Schrodinger equations, which arc not 
conservation laws. In particular, the pseudo-conformal identity and Morawetz' identity are im- 
portant in scattering theory for repulsive <i>'s. Another identity referred to as variance identity is 
an ordinary differential equation. It can be applied to exhibit blow-up of solutions of ljIIE|) for 
non-subcritical, attractive two-body potentials. For example, choosing A = and i'^{x) — — |a;|^°" 
the variance identity becomes 

\x\'\m\'dx^l6nm + {4-2a) f {Ix^ * m)\^)m)\' , (3.16) 

Clr JR3 Jr3 

provided that tpo G {?/' G X : \x\'^\'ip\'^ G i-'^(M'^)}. Since the left hand side is the second derivative 
of a nonnegative function which cannot be strictly concave for all t > 0, one can infer that ^(t) 
can only exist over a finite time interval if the right hand side is bounded from above by a negative 
constant. This happens for tr > 2 (i.e., <& is non-subcritical) and initial data with 7i[V'o] < 0- 
We refer to |Uaz03j for a detailed exposition of such identities and their applications to nonlinear 
Schrodinger equations. 



4 Hartree Solitons 

In this section we discuss solitary wave solutions for (jHEp . which we call Hartree solitons. By 
using variational methods, we derive existence and orbital stability of Hartree ground state soli- 
tons, for attractive, subcritical <I>'s. See Theorems [3 and |3| below. Orbital stability for local, 
subcritical nonhnear Schrodinger equations has only been established for special classes, e.g., for 
pure power nonlinearities and some variations thereof. We also demonstrate how linear stability 
can be understood from the variational proof of stability. Finally, we touch upon the issue of 
asymptotic stability, which involves a detailed study of dispersive estimates of the linearization 
around a Hartree soliton. Some open problems are proposed. 

Plugging the ansatz ^(f, x) = e"^*Q(a;), with e M, into (|HE|I leads to the nonlinear eigenvalue 
problem: 

-l^Q + \VQ + y{^*\Q\^)Q^~LoQ, I 
Q e X, Af[Q] . 

Recall that X denotes the energy space introduced in H2.13|l . We exclude the trivial solution, 
Q = 0, by requiring that N[Q] = ||Q||i = iV > 0. 
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We note that there are two ways of treating problem (|4.1(l : Either we fix N and try to solve 
for Lo = lon and Q — Qn, or we fix lj and try to solve for N — N^^ and Q = Q^. In physical 
applications it is sometimes convenient to rephrase the first way by fixing N as an overall constant 
and then look for solutions Q — Q^, and lo — lOi, depending on the coupling constant v. 

If there exists a solution Q — Q^^ of 1)4. 1|1 . with w G / C M, then 

:=e''"*Q^(x), uo € I , (4.2) 

is a family of solitary wave solutions for l)IIE() . which we refer to as Hartree solitons. Provided 
that A = holds, the family of Hartree solitons comprises the functions 

^,{t,x) e*(^-(--'*)-^'*+^+"*)Q^(a: - vt - d) , 
a (d, V, 7, tj) e X X [0, 27r) X / , 

obtained by applying Galilei transformations to ip^^ (t, x) ; spatial rotations can be omitted if Q is 
a radial function. See Theorem |21 below for sufficient conditions for Q to be radial. 



4.1 Existence of Hartree Ground States 

We relate the nonlinear eigenvalue problem 1)4.1(1 to the following variational problem 

^(A^) := inf{7^[V'] : V e Sjv} , (4.4) 

where 

Sat = {iA e X : A/'[^/;] = TV} . (4.5) 

If a minimizer exists then it is a solution of ()4.1|l . provided that 7i is C^. Such energy minimizing 
solutions are called Hartree ground states. To arrive at a general statement on the existence 
of Hartree ground states, we require the following definition: We say that $ is subcritical and 
vanishing at infinity if 

$ e LP' (R3) + (]g3) ^ fQj, gQj^g 3/2 < pi < p2 < oo , (4.6) 

and likewise for V. 

Theorem 2. [Existence of Hartree Ground States] Without loss of generality we set v — 1 
in \4-l\j - Suppose that $ is real-valued, radial, subcritical and vanishing at infinity, with $ < 0, 
$ ^ 0. Let V be real-valued and assume that one of the following assumptions is satisfied. 

a) A = (i. e., no external potential) . 

b) V is subcritical and vanishing at infinity with V < Q on a set of positive measure (i. e., V is 
a "localized trap"), A = 1. 

c) V L^^^iM.^), V > 0, and V{x) — > oo as |x| —> oo (i. e., V is a "confining trap"), A = 1. 

Then there exists a constant N^, > such that, for every N > N^,, problem j^. j| j has a Hartree 
ground state solution, Q — Qn, which can be chosen to be strictly positive. 

In addition, if is nondecreasing then Q is radial (about some point) and nonincreasing, 
provided that either a) holds, or b) holds with V radial and nondecreasing. 
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Remarks. 1) Under condition a) one can show that iV* = 0, for <i>(x) = — |a;|^'^, with < a < 2. 
However, for two-body potentials of short range (e. g., of van der Waals-type) N^, is strictly positive. 
With regard to condition c), it is easy to see that = 0, for confining Vs. 

2) By applying standard methods we can derive decay and regularity properties for Qn de- 
pending on V. If, for instance, condition a) holds we can infer that Qn G C°°(M.^), thanks to 
elliptic bootstrapping, and that Qn exhibits exponential decay at infinity. 

3) Theorem 121 can be rephrased as follows: For iV = 1 there exists a critical coupling, > 0, 
such that, for all i' > v^, there exists a Hartrcc ground state Q — Q^. 

Sketch of Proof. The assumptions of Theorem El imply that, for every e G (0, 1), there exists a 
constant C^^n > such that 

EiN)>^\\\7m-C,,N. (4.7) 

In particular, this shows that E{N) > — oo holds for all N < oo. To prove existence of Hartree 
ground states, for N sufficiently large, we apply the direct method of the calculus of variations. 
That is, we consider minimizing sequences, (V'n)n>Oj in Sat with 

n[tpn] E{N) , as 71 oo , 

and try to find a strongly convergent subsequence in X. Depending on conditions a), b), or c) in 
Theorem 12 we can show relative compactness of every minimizing sequence as follows. 

Relative Compactness under Condition a). Due to the assumptions on $, we can apply the 
concentration-compactness method developed by P.-L. Lions |Lio84al ILio84bj for translation- 
invariant functionals. In these references, it is shown that the strict sub-additivity condition 

E{N) < E{a) + E{N - a) , for all a G (0, N) , (4.8) 

is sufficient and necessary to show that every minimizing sequence, ('0„)n>o, of 7i['0]|^gSjv is 
relatively compact in X = H^(M.^) up to translations. Thus, there exists a sequence, {yn)n>o, 
in such that {tpn{- + yn))n>o has a subsequence that is strongly convergent in H^{M.^). In 
particular, a minimizer exists if this property holds. 

We must show that 1)4 .8(1 holds for all N > N^, for some constant iV, > 0. This is accomplished 
as follows. First, it is easy to see that 

E{N) < , for all TV > . (4.9) 

Moreover, we infer from the scaling behavior 

E{^N)^m\^ f [|V7/^|2 + ^($*|^|2)|V;|2]dx: VeS^I (4.10) 

that there exists some TV, > such that 

E{N) < , for all N > , (4.11) 

and that 

E{m) < ^E{N) , for aU N > and > 1 . (4.12) 

The strict sub-additivity condition (|4.8|) for N > N^, then follows by elementary arguments. 

Relative Compactness under Condition b). To show that every minimizing sequence for N > 
has a strongly convergent subsequence in X = H^{M.'^) (up to translations), we can proceed 
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similarly as in a). See jLio84al |Lio84b| for a modified sub-additivity condition for translation- 
dependent functionals. 

Relative Compactness under Condition c). If is confining then it easy to see that, for every 
e > 0, there exist constants Re,n^ > such that 

/ |^„Pd.<s, forall n>n.. 

J\x\>Rc 

This fact allows us to pass to a strongly convergent subsequence for every minimizing sequence in 
X. 

Strict Positivity. Note that "^[IV"!] ^ "'^[V'Ij because ||V|'i/'|||2 < ||Vi/'||2- Thus, we can choose 
Qn to be real- valued and nonncgative. Introducing the Schrodinger operator 

H:=-A + W, where W -.^ V + {<P * {Qn^) , (4.13) 

we have that Q^r is an eigenfunction, i. e., 

HQn^-ujnQn, (4.14) 

which can be written as follows 

e-^Q^r ^e-^^QAT. (4.15) 

The assumptions of Theorem El imply that W+ :— max{T4^, 0} G ^[^^.(M'^) and that :— 
min{iy, 0} is subcritical. With results from |FS75j we find that e^^ is positivity improving, i.e., 
/ > 0, with / ^ 0, implies that e^^ f > 0. Hence Qn{x) > almost everywhere. 

Radial Symmetry. The statement about radial symmetry of Qj\j in Theorem [21 follows from 

'H[Q*N]<n[QN], (4.16) 

provided that $ = $* and V — V*, where /* denotes the symmetric rearrangement of a measur- 
able function / that vanishes at infinity; see, e.g., jLLDlj for rearrangement inequalities. □ 

Remarks. 1) In physical terms, condition 1)4. 8|) says that Hartree ground states are stable against 
fission: A ground state configuration having mass N has an energy strictly smaller than the sum 
of energies of two configurations with mass a and N — a, respectively. Hence there is always some 
strictly negative binding energy. In mathematical terms, condition (|4.8|1 prevents every minimizing 
sequence from either converging to zero locally everywhere, or splitting up into lumps receding 
from each other. 

2) If N is sufficiently large then nonlinear excited states (which are not strictly positive) exist, 
too. They are critical points of Hlip] restricted to Sat; see |Lio80| for a critical point theory. 

Theorem El does not answer the question of uniqueness of Hartree ground states, i. e., whether 
all minimizers of 7i[V'] on S^r are just the orbit of some Qn under the group of symmetries leaving 
7i[tp] invariant. By using a remarkable trick, E. Lieb |Lie77| has proven uniqueness of Hartree 
minimizers for H2.16|l . meaning that all minimizers are equivalent up to phases and translations. 
But the general case is not understood, yet. Thus, we propose. 

Problem 3. Prove uniqueness of Hartree ground states on Sn, up to symmetries o/7i['0], for a 
general class of $ 's and A = 0. 
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Remarks. 1) Uniqueness of positive ground state solutions for nonlinear Schrodinger equations on 
R*^ with local nonlinearities of the form is a known fact, for < a < 4/((i — 2); see |Kwo89| . 

The proof given there relies, however, on Sturm comparison theorems and ordinary differential 
equation techniques that cannot be directly applied to nonlocal equations such as (|4.1() . 

2) If an external potential is present (i. e., A 7^ 0), then one can also study perturbative regimes 
to show existence and uniqueness of Hartree ground states. Assuming that —A + XV has a nor- 
malizable ground state, one can prove uniqueness of Hartree ground states for ?i['0]|s„, with N 
fixed and v sufficiently small; see, e.g., |AFG+02] for such results. 

Finally, we wish to draw attention to the phenomenon of symmetry breaking: If, for example, V 
is a double- well potential then, for suitable $ and sufficiently large z/, two distinct Hartree ground 
states exist; i.e., uniqueness is lost |AFG"'"02] . (Of course, translation invariance is always broken 
when A = 0, as soon as a non-trivial Hartree ground state exists.) 

4.2 Orbital Stability of Hartree Ground State Solitons 

Next, we study orbital stability of Hartree solitons arising from the Hartree ground states con- 
structed in Theorem 13 For this purpose, we introduce the set of Hartree ground states, Mn, for 
N > N,, by 

Mn ■■= {Q (^Sn -.niQ] ^ mm H[^P]} , (4.17) 
and define the distance function 

d(V'):= inf U-QWx- (4.18) 
QeMN 

We note that if we have uniqueness of Hartree ground states (up to phases and translations) , then 

d{^j) = inf inf - e'^Q(- - y)||x (4.19) 

for a fixed Q £ Mn- Recall that uniqueness of Hartree ground states is known for equation (|2.16(l . 
which models a non-relativistic boson star; see jLie77j . 

We say that A^^f is orbitally stable in X (energy space) if, for every e > 0, there exists S > 
such that, for all initial data ipQ d X, 

d{ipo)<S implies that d{ip{t)) < e , for all t>0, (4.20) 

where ip{t) denotes the solution of (|HEp . With the proof of Theorem|2]at hand, it is easy to prove 
the following result. 

Theorem 3. [Orbital Stability of Hartree Ground State Solitons] Assume that the hy- 
potheses of Theorem\^are satisfied (and if V is confining that condition b) of Theorem^ holds). 
Then orbital stability of the set of Hartree ground states, M. n , holds for all N > N^. 

Remark. If the external potential vanishes (i.e., A = 0), then Theorem |31 also guarantees orbital 
stability of "boosted" Hartree ground state solitons; see equation H4.3|l . This fact follows by using 
the Galilei invariance of l)HE|l . 

Sketch of Proof. The proof of Theorem O makes use of conservation laws and the relative compact- 
ness of every minimizing sequence; the idea goes back to |(]L82| . It proceeds as follows. Recall 
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from Theorem n that (|HE|I is globally well-posed in X . To simplify the argument, we restrict our 
attention to initial data with A/'['0o] = N (the general case will then follow). Let us assume that 
orbital stability does not hold. Then there would exist an e > and a sequence (u„)n>o in Sat 
with 

d(un) ^ , as n ^ oo , (4-21) 
such that, for the solution tpnit) of (|HE|I with initial datum u„, 

d{ipn{tn)) > £ , for aU n > (4.22) 

and a suitable sequence (i„)„>o of times. Let us define w„ :— ijjn{tn)- Since charge and energy are 
conserved, i. e., 

M[vn] = AA[u„] , n[vn] = n[u„] , (4.23) 

the sequence (w„)„>o is also a minimizing sequence for 7i[V']|sjv- have already seen in the 

proof of Theorem 121 every such minimizing sequence contains a strongly convergent subsequence 
in X (up to possible translations). This contradicts (|4.22|) and proves Theorem O □ 



Orbital Stability via Linearization 

There is another way to discuss orbital stability involving the linearization around a Hartree soli- 
ton. In fact, a rigorous analysis, initiated by M. Weinstein |Wei85| for local nonlinear Schrodinger 
equations and extended by M. Grillakis, J. Shatah, and W. Strauss jCSSQOj . shows that the 
stability condition 

^AA[Q^] > (4.24) 

is sufficient for orbital stability of Hartree ground states. Recall that Q = denotes the solution 
of H4.1() . for a given uj, choosing N = Af[Quj]. Condition H4.24II is popular in the physics literature 
on solitons. 

Under the assumptions of Theorem [21 we now demonstrate why condition 14.24(1 implies the 
strict sub-additivity condition (|4.8|l and thus orbital stability. It turns out that the relation 
between these two criteria is essentially the same as the one between the inequality /" < and 
the strict concavity of a function /. 

To convince ourselves of this claim, we assume that Q = Qn and Q = can be differentiated 
with respect to A'^ or a; up to second order. We observe that the nonlinear eigenvalue problem 
(|4.1|) with constraint Af[Q] = N is just the Lagrange equation with multiplier —uj/2, i. e., 



■H'[Qn]^-^^'[Qn]- (4.25) 



Thus, for all iV > iV*, we find that 



(4.26) 



_ lu(W[Qn] _ ujdN _ UJ 
~ ^2 dN ^ ^2 diV ^ ^2 ' 

The stability condition ((4.24() then implies that 

i?"W = -^^<0. (4.27) 

Hence E{N) is a strictly concave function for all N > N^,. The strict sub-additivity condition 
()4.8|l is a direct consequence of this fact. 
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4.3 Instability and Blow-Up 

Using the concentration-compactness method to prove orbital stabihty is quite elegant. It is, 
however, not suited to disprove orbital stability if $ is not subcritical, because it is possible that 
E{N) — —oo, i.e., that the Hartree energy functional is not bounded from below on S^v (hence 
A^AT is empty). If this happens, linearization around Hartree solitons, defined to be critical points 
o{Ti.[ip] on Sat, appears to be indispensable. In fact, it can be shown that the instability condition 

< (4-28) 

is sufficient for orbital instability for Q = [GSSQOj . This means that there exists an e > 0, but 
no (5 > such that d{ipo) < S implies that d{%l){t)) < e, for all t > 0, where d{-) is defined as 

:= inf inf - e'^Q{- ~y)\\x- (4.29) 
As an example, we illustrate instability of Hartree solitons for the non-subcritical case, 

i^tVj = -AV-- (la;!"" * IV'HV', with 2<(T<3, (4.30) 
where non-trivial solutions, Q = Q^^, of the equation 

-AQ- {\x\-^ ^\Q\^)Q = -LuQ, (4.31) 
are known to exist; see |Lio80| . A calculation shows that 

M[Q^] = Cuj^~^ (4.32) 

for some positive constant C — C{a). Hence, the instability condition H4.28|l holds for 2 < ct < 3. 
That orbital instability also occurs for a = 2 (where condition H4.28|l fails) can be demonstrated 
with the help of the variance identity (|3.16|) as follows. Multiplying equation H4.31|l by Q and 
integrating over space, and doing the same with x ■ VQ, leads to the identity 

n[Q] = 0. (4.33) 

Choosing an initial datum ipo := (1 -I- e)Q, we find that 7i[V'o] < 0, for all £ > 0. According to 
the variance identity (|3.16|) and the discussion following it, the corresponding solution of (|HE|I 
exhibits blow-up in finite time; (note that ipQ has finite variance because of its exponential decay). 
Choosing e sufficiently small we can approximate Q by -00 arbitrarily well in X = iJ^(R'^). This 
shows instability of ^(t,x) = e'-'^^Q{x) in (|4.3()() . for a = 2. 

4.4 Asymptotic Stability 

Refining the stability analysis leads one to the issue of asymptotic stability of Hartree solitons. 
One asks whether in some suitable norm, || • 

||V'(i,-)-e*'''*'"^Q^(0(--yW)||y ^0, as < ^ +oo , (4.34) 

for all initial data tpQ sufficiently close to a soliton configuration and some functions •d{t^x)^Lj{t)^ 
and y{t) with the properties that 

i^(^) a;) iwoo • (a; - j/oo) - 71^^^ + 7cx) + Wooi , iu[t)^uJoo, y{t) ^ Voo + v^t , (4.35) 
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as t — > +00. A generalization of the notion of asymptotic stability to multi-soliton configurations 
can be formulated easily. 

We briefly outline the history of work on asymptotic stability for nonlinear Schrodinger equa- 
tions. First results date back to the work of A. Soffer and M. Weinstein |S W90j for local nonlinear 
Schrodinger equations with suitable external potentials and further restrictive conditions. Their 
work has been motivated, in part, by conjectures of J. F. and T. Spencer. (For nonlinear Klein- 
Gordon equations with spatial inhomogeneities, first results appeared in |Sig93| and were extended 
in |SW99| ) . The first work on asymptotic stability for translation- invariant nonlinear Schrodinger 
equations (also under very restrictive conditions) was presented by V. Buslaev and G. Perel'man 
rRP92) . Both works were significantly extended in jSW92l ISWHrI rRP95l Th^di^ ITYn2bl rTYn2al 
ITY02cl IGNTM) . See |PW97I IWedOO) for a center-manifold approach. Recent results concern- 
ing scattering of multi-soliton configurations for local nonlinear Schrodinger equations have been 
established in |Per03[ 1^5^03) : (see also |FTY02j for partial results). 

As far as extending such methods and results to Hartree-type nonlinearities is concerned, one 
expects that there are natural conditions on $ allowing for asymptotic stability. Such results 
would have interesting physical applications. 

To propose a problem concerning asymptotic stability of Hartree solitons, we assume that 
A = and consider the linearization of (|HEp around a Hartree ground state Q^^. Making the 
ansatz 

i;{t, x) = e'^'lQM) + h{t, x)] (4.36) 

leads to the linear equation 

idth^L^h, (4.37) 

where 

L^h ^ ^Ah + Luh + {<P * \Q^f)h + 2Q^{^ * {Q^{h + h))) . (4.38) 

Note that is not C-linear and nonlocal. As an immediate consequence of the symmetries of 
l|HEp with respect to U{1) gauge-transformations and spatial translations one finds that 

L^iQ^^O, L^djQ^^O, where j = l,2,3. (4.39) 

This observation leads to the following problem. 

Problem 4. Find a class of two-body potentials $ with the following properties. 

a) One has dispersive estimates for L^. 

b ) The null space of is given by 

Null (L^) = spauR {iQcj, djQ^} , where j ^ 1,2,3, 

as expected on the basis of the symmetries of I^HF]) . 

Remark. Note that the local potential, $ * IQc^P, in exhibits power-law decay if <i> has power- 
law decay. In contrast, local nonlinearities lead to exponential decay, thanks to the exponential 
decay of Q^^. Exponential decay of the background potential in the linearization around a soliton 
plays a central role in the derivation of dispersive estimates in |RSS03| . 



20 



5 Point-Particle Limit 



In this section, we discuss the point-particle (Newtonian) limit of the Hartree equation (jHE(l and 
show how Newton's equations of motion for a system of point-particles emerge. 

5.1 Basic Ideas and Result 

To exhibit the point-particle limit, we consider the following situation. The external potential is 
given by 

XV{x) := V'^'\x) := W{ex) , (5.1) 

where W{x) is a smooth, bounded function and e > denotes a parameter. Furthermore, we 
assume the two-body potential to be given by 

^(x) = <i>six) + ^eiex). (5.2) 

Here $s is a radial, short-range, attractive potential, whereas ^£ is smooth, radial and may be 
either attractive or repulsive. In addition, $^ is allowed to be of long range, e.g., the Coulomb 
potential. 

As a starting point, let us consider a configuration of Hartree ground states for $ = $s (i. e., 
e = 0). We pick k numbers A^i, . . . , Nk > N^(^s) and Hartree ground states, Qnj, minimizing 
Ti-lip] on Snj, for A = and <i>^ = (see Theorem |2l. We consider an initial configuration of k 
solitons 

{Qiv,(x-<Z,)}'^^ (5.3) 
centered at points 9j G K'^, 1 < j < k. These solitons have a characteristic diameter given by 

Lsol := max ■ /— . (5.4) 

j = l,...,k y LUj 

The external potential, V'-'^^{x) = W{ex), defines another characteristic length 

esuplVW^I 

Choosing e > sufficiently small, such that 

Lsol sup |VW| 



<1, (5.6) 

one finds that the external potential varies little over the typical spatial extension of a soliton, 
i.e., the solitons are almost "point-like". Finally, we adjust the initial positions of the solitons to 
be far separated from each other, i. e., 

e>—^ r. (5.7) 

mm^^j \qi - qj] 



Having specified our initial configuration, we now look for solutions of (jHEjl of the form 

k 

{t, x)^J2 (*) - * + '^^'^ ^) (5.8) 
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obeying, for all times t with |t| < 0(e the following equations 

0j {t, x) = mqj {t) ■ [x - qj (t)] + ?9j (t) , (5.9) 

\N,{t)\ = o{e) , (5.10) 

\\h^'\t,-)\\H^ ^oie), (5.11) 

where we reinstall the boson mass m. A formal calculation with the action 13.8|l . using the ansatz 
(|5.8|l . yields 

"J=i (5.12) 



EW.-<f^ (£(g.(t)-g,(i)))+o(£) 



i<j 



for r = ©(e"^). For the action to be stationary at -0 = V'^^^ meaning that ip'^'^^ solves (|HE|) . 
the set of trajectories, {(/j(i)}j=iJ must apparently satisfy Newton's equations of motion for k 
point-particles with mutual interactions given by $f and subject to an external potential W, up 
to small friction forces: 

mqj{t) = -eiVW)ieq,{t)) - e E W - + W , (5.13) 

for = 1, . . . . fc and times t with |i| < 0{e^^). The inhomogeneities, aj{t), obey the estimate 

\a,{t)\^o{e), (5.14) 

and are interpreted as "radiation damping" or "friction" . 

All these considerations can be formulated as a theorem for fc = 1 jFTY02| . but the methods 
can be extended to general k. A more systematic approach to such problems can be found in 
|FGJS04| . where use of the Hamiltonian formalism and the associated symplectic structure is 
made. 

6 More Open Problems, Conclusions, and Outlook 

Recall the open problems already mentioned above: Problem ^ (mean- field spectral analysis), 
Problem|21 (relaxation into coherent states), Problem|31 (uniqueness of Hartree ground states), and 
Problem^ (dispersive estimates for asymptotic stability). In addition to these problems, we wish 
to draw the reader's attention to the following open problems, which we think are worth studying. 

Integrable Motions 

Consider ljHE() with an external harmonic trap: 

i9tV' = -AV'-l-|2:|V+(**lV'nV'- (6.1) 
One easily verifies that if Q{t,x) solves (|6.1f) then another solution is given by 

^Pit, x) := e^[f=W-"-^*lg(i, X - x,{t)) , (6.2) 
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where {pc{t), Xc{t)} solve the classical equations of motion of a harmonic oscillator in R'^, i. e., 

Mt)^-2x^{t), Xcit) = 2p^{t) , where E := {p^{t)f + {x^{t))^ ^ const. (6.3) 



Thus, we can construct solutions of 1)6.1(1 corresponding to Hartree solitons "moving" along orbits 
of an integrable point-particle system. 

The following generalization is of interest; Consider an external potential V that is radially 
symmetric. Then the Hamilton function 



describes an integrable Hamiltonian system. Choose V such that the system satisfies strong 
anisochronicity conditions. Using "singular" KAM theory, attempt to construct integrable motions 
of Hartree solitons, for iV sufficiently large, whose centers of mass follow orbits of (|6.4|l . 

Formation of Hartree Solitons 

In view of physical applications, e.g., to boson stars, it would be interesting to gain a deeper 
understanding of the formation of Hartree solitons. Consider, for example, a configuration of 
two far-separated solitons with masses Ni,N2, initial positions 91,(72, and initial velocities vi = 
91(0), ^2 = 9(0), as described in equation 1(5. 8|l . with h^'^\x,0) = 0. Assuming that the two-body 
potential, $f , is purely attractive and of short range, we choose {Ni, qi,Vi}'^^i in such a way that 



i.e., we consider a bound configuration of two solitons. What happens to this system as t ^ oo7 
Intuitively, we expect collapse into a single soliton moving inertially through space (for A = 0), 
perturbed by some outgoing dispersive radiation tending to zero at the free dispersion rate. There 
is numerical evidence for such a behavior, but analytical results are still lacking. 

We also mention recent results on soliton formation derived in |Tao04j . The ideas developed 
there appear applicable to l)HE() with attractive two-body potentials of short range. 

Random External Potentials 

In the study of transport phenomena, it is interesting to consider (IHE|I with a random external 
potential, V, whose distribution is homogeneous and has correlations of short range; see |AF88I 
IAFS88I IBW04| for local nonlinear Schrodinger equations with random external potentials, and 
|FSW86j for KAM methods applicable to some random systems with infinitely many degrees of 
freedom. 

Renormalization Group Methods 

Recent work on gravitational collapse incorporates ideas from the theory of critical phenomena 
and renormalization group methods; see, e.g., |Gun99| . It would be of interest to apply such 
methods to the semi-relativistic Hartree equations (I2.17|l modeling a boson star; see |Kup95| for 
a general discussion of renormalization group methods in the analysis of blow-up of solutions of 
nonlinear partial differential equations. 



/ + Vix) 



(6.4) 



Nivl + N2VI + NiN2^e{e{qi - (72)) < , 



(6.5) 
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